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Our Local Group of galaxies is moving with respect to the 
cosmic microwave background (CMB) with a velocity' of 
Vems = 631+ 20kms— and participates in a bulk flow that 
extends out to distances of -20,000 km s~’ or more?“, There 
has been an implicit assumption that overabundances of gal- 
axies induce the Local Group motion®”. Yet underdense regions 
push as much as overdensities attract®, but they are deficient 
in light and consequently difficult to chart. It was suggested a 
decade ago that an underdensity in the northern hemisphere 
roughly 15,000 km s~ away contributes significantly to the 
observed flow’. We show here that repulsion from an under- 
density is important and that the dominant influences causing 
the observed flow are a single attractor — associated with the 
Shapley concentration — and a single previously unidentified 
repeller, which contribute roughly equally to the CMB dipole. 
The bulk flow is closely anti-aligned with the repeller out to 
16,000 + 4,500 km s~’. This ‘dipole repeller' is predicted to be 
associated with a void in the distribution of galaxies. 

The large-scale structure of the Universe is encoded in the flow 
field of galaxies. A detailed analysis of the flow uncovers the rich 
structure manifested by the distribution of galaxies, such as the 
prominent nearby clusters'*"", the Laniakea supercluster’* and the 
Arrowhead mini-supercluster'®. A close correspondence between 
the observed density field, derived from redshift surveys, and the 
reconstructed three-dimensional (3D) flow field has been estab- 
lished out to beyond 100 megaparsecs and down to a resolution of 
a few megaparsecs'®. Yet the flow contains more information on 
distant structures, as found from tides and from continuity across 
the zone obscured by the Galactic disk, the ‘zone of avoidance’!’. 
The Cosmicflows-2 dataset of galaxy distances’® provides reason- 
ably dense coverage to R% 10,000 km s™ (distances are expressed in 
terms of their equivalent Hubble velocity). 

The linear 3D velocity field is reconstructed here from the 
Cosmicflows-2 data by the Bayesian methodology of the Wiener 
filter and constrained realizations (see Methods). The Wiener filter 
is a Bayesian estimator that assumes a prior model: here it is the 
ACDM model. It is a conservative estimator which balances between 
the data and its errors and the assumed prior model. Where the data 
are weak, the Wiener-filter estimation tends to the null hypothesis 
of a homogenous universe: that is, a vanishing peculiar velocity. Yet 
it uncovers structures out to ~16,000kms™! at the extremity of the 
Cosmicflows-2 coverage. The variance around the mean Wiener- 
filter estimator is sampled by the constrained realizations. 

The Wiener-filter flow field is used to construct the cosmic web, 
defined here by means of the velocity field” (V-web; see Methods) 
This is done by evaluating the velocity shear tensor on a grid and 
counting the number of its eigenvalues above a threshold value: 3 
corresponds to a knot, 2 to a filament, 1 to a sheet and 0 to a void. 
In the linear regime, the flow is proportional to the negative of the 


gravitational field; hence it constitutes a gradient of a scalar poten- 
tial. Figure 1 shows the large-scale structure out to a distance of 
16,000 kms! in a plane that contains the Local Group, the Shapley 
attractor and the dipole repeller. Three different aspects of the flow 
are depicted: streamlines, the V-web and the velocity potential. 

When describing the gravitational dynamics in co-moving coor- 
dinates, by which the expansion of the Universe is factored out, 
underdensities repel and overdensities attract. The velocity field is 
represented here by means of streamlines (see Methods), the sources 
and sinks of which are the attractors and repellers of the large-scale 
structure (see Fig. 1 and the Supplementary Video). Figure 2 shows 
a 3D visualization of the streamlines in a box of length 40,000 kms! 
centred on the Local Group. The streamlines of the left plot of Fig. 2 
either converge onto an attractor located roughly at [— 12,300, 7,400, 
—300] kms“ or cross out of the box. The plot uncovers the existence 
of a repeller at the upper right-hand side of the box — a region from 
which streamlines diverge. Repellers are best manifested by the anti- 
flow, namely the negative of the velocity field. The right plot of Fig. 2 
depicts the convergence of the anti-flow streamlines onto a repeller 
at supergalactic coordinates [11,000, —6,000, 10,000] kms~'. The 
Wiener-filter reconstruction detects a single attractor and a single 
repeller, the Shapley attractor and the dipole repeller. (The accom- 
panying Supplementary Video provides a further visualization.) 

Multipole expansion provides a different insight into the nature 
of the observed flow. The dipole (bulk flow) and quadrupole (shear 
tensor) moments are evaluated here by spherical top-hat window 
functions'''*-”° centred on the Local Group with variable radii, R. 
The telltale signature of a single dominant attractor or repeller is 
the close alignment of the dipole with the expansion eigenvector 
of the shear tensor and a degeneracy of the other two eigenvalues/ 
eigenvectors. The observed flow is clearly not dominated by either 
a single attractor or repeller. In the following, we emphasize the 
directional aspects of the dipole and shear eigenvectors, which are 
robustly recovered by the Wiener filter. Figure 3 presents an Aitoff 
projection of the following directions: (1) the dipole repeller; (2) the 
Shapley attractor; (3) the CMB dipole and its anti-apex; (4) the bulk 
velocity of top-hat spheres of R = (2,000, 3,000,..., 15,000) kms“, 
V ur(R), of the Wiener-filter reconstructed flow field; (5) the three 
eigenvectors of the shear tensor (ê; i= 1,2,3) of the Wiener filter 
field. The figure shows the strong anti-alignment of the bulk flow 
out to 15,000kms~! with the dipole repeller. Beyond that radius, 
the bulk flow loses its coherence, as the scatter in direction steadily 
increases. The eigenvector of the shear tensor that reflects the direc- 
tion of maximal expansion (€;) is aligned with the direction of the 
Shapley attractor out to R = 7,000km s™. Supplementary Fig. 2 
further presents the mean and the scatter around the cosine of the 
angles formed between the bulk velocity and the dipole repeller, 
Houlk(R) = cos(Vbuk Rar), and between e, and the Shapley attractor, 
He;(R) = cos(é3(R), R shapley)- 
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Figure 1 | A face-on view of a slice 6,000 kms“ thick, normal to the direction of the pointing vector f = (0.604,0.720, — 0.342). Three different elements 
of the flow are presented: mapping of the velocity field is shown by means of streamlines (seeded randomly in the slice); red and grey surfaces present the 
knots and filaments of the V-web, respectively; and equi-gravitational potential (@) surfaces are shown in green and yellow. The potential surfaces enclose 
the dipole repeller (in yellow) and the Shapley attractor (in green) that dominate the flow. The yellow arrow originates at our position and indicates the 
direction of the CMB dipole (galactic longitude |= 276°, galactic latitude b = 30°). The distance scale is given in units of kms“. 


Two aspects of the bulk velocity — the anti-alignment itself 
(Mpur(R) = —0.96 + 0.04) and its distance scale (R = 16,000 kms“) 
— strongly corroborate the Wiener-filter finding of the dipole 
repeller and its dominant role in dictating the observed flow. 
The direction and distance coincide with the position of the dipole 
repeller — a non-trivial occurrence. It is interesting to study the 
shear tensor. The expansion eigenvector is closely aligned with 
the direction to the Shapley attractor out to Ræ 7,000kms™', the 
distance of the foreground (Norma-Centaurus—Hydra) Great 
Attractor, located at the bottom of the Laniakea basin of attrac- 
tion at supergalactic coordinates (—4,700, 1,300, 500) km s™!. It is 
the combined mass distribution within the Laniakea and Shapley 
superclusters that dominates the tidal field, with the inverse cubic 
distance dependence of the tidal interaction tipping the balance at 
our location towards the Laniakea/Great Attractor. 

Our main findings are tested against statistical and systematic 
uncertainties. There is no doubt about the existence of the Shapley 
concentration, and therefore we focus our attention mostly on the 
dipole repeller. The strong support for the existence of the dipole 
repeller comes not only from its close alignment of the bulk velocity 
but also from the small scatter around the mean Wiener-filter value, 
Houl R) =—0.96 + 0.04 for R x 16,000 km s™ (Supplementary Fig. 2). 
Assuming that the dipole repeller is the dominant structure that 
determines the direction of the bulk flow, the scatter in ,,,.(R) 
can be translated to uncertainty in the position of the repeller, 
ARpr Z 4,500 kms"! (see Methods). The basins of repulsion and 
attraction around the dipole repeller and Shapley attractor, out to 
a distance of ~8,000 km s~’, contribute fairly evenly to the velocity 
of the Local Group. At 8,000 km s™', the repeller and the attractor 
contribute 59 + 26 and 67 + 27kms", respectively, to the CMB 
dipole (see Methods). Next, the robustness of the dipole repeller is 


tested against subsampling of the data. These consist of cuts either 
by distance (6,000, 8,000 and 10,000 km s~’), or by galaxy and data 
type (see Methods), corresponding to a degradation in the quality 
of the data by volume coverage, number of data points and magni- 
tude of errors (see Methods). All subsamples considered here locate 
the dipole repeller in an underdense region and recover a basin 
of repulsion that pushes the Local Group in the direction of the 
CMB dipole. 

The general picture that emerges here is of a complex flow that 
cannot be explained by a simple toy model, yet the main structures 
that shape the observed flow can be identified. The Wiener filter 
recovers a flow dominated by a single attractor and a single repeller, 
which roughly equally contribute to the CMB dipole. The role 
played by the Shapley attractor is not surprising; the earlier find- 
ings on influences beyond the Great Attractor®”*?!” suggested 
it. The existence of the dipole repeller was only vaguely hinted at 
before. A study of the all-sky distribution of X-ray-selected clus- 
ters uncovered a significant underdensity of clusters in the north- 
ern hemisphere roughly 15,000 kms~! away’. It suggested that this 
underdensity may be as significant as the overdensity of clusters in 
the southern hemisphere in inducing the local flow. Earlier exami- 
nations of galaxy peculiar velocities found a north-south anisot- 
ropy in (galactic) y-component of the velocities’ and found” that 
the sources responsible for the bulk flow are at an effective distance 
>30,000 km s~. Here, the source of the repulsion is identified for 
the first time. The dual dominance of the dipole repeller and the 
Shapley attractor is the main new finding of this study. The strong 
anti-alignment of the CMB dipole with the dipole repeller out to a 
distance of 16,000+4,500 kms~! suggests the possible dominance 
of the repeller over the attractor. The predicted position of the 
dipole repeller is in a region that is as yet poorly covered by existing 
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Figure 2 | A 3D view of the velocity field. It is shown here by means of the flow streamlines (in black-blue, left panel) and of the anti-flow (in yellow-red, 
right panel). Anti-flow is defined here by the negative (namely, the reverse) of the velocity field. The same streamlines are seeded on a regular grid 

and are coloured according to the magnitude of the velocity. The flow streamlines diverge from the repeller and converge on the attractor. For the 
anti-flow, the divergence and convergence switch roles: they diverge from the attractor and converge on the repeller. The knots and filaments of the 
V-web are shown for reference. Cartesian supergalactic coordinates (SGX, SGY, SGZ) are assumed here. (For a 3D view, look at the accompanying 


Supplementary Video, at time 00:56-01:28.) 


redshift surveys. We predict the dipole repeller to be associated with 
a void in the distribution of galaxies. 

In the linear regime of gravitational instability, repellers are as 
abundant and dominant as attractors. Yet, observationally, repel- 
lers are much harder to identify than attractors. The association of 
repellers with underdensities renders them strongly deficient in gal- 
axies, in general, and clusters of galaxies, in particular, and thereby 
makes their direct detection challenging. Our use of peculiar veloci- 
ties as tracers of the large-scale structure overcomes that observa- 
tional hindrance and unveils the existence of the new extended 
structure located at a=22h 25min, 6 = +37° (galactic longitude 
1=93°, galactic latitude b=—17°; supergalactic latitude and longi- 
tude SGL = 332°, SGB = 39°) that we call the ‘dipole repeller’ 


Methods 


Cosmicflows-2 dataset. ‘The present study is based on the second release 
catalogue” of galaxy distances and peculiar velocities, Cosmicflows-2, which 
extends sparsely to recession velocities of 30,000 km s~ (redshift z% 0.1). It consists 
of 8,161 entries with high density of coverage inside 10,000kms~'. Here we used 

a grouped version of the Cosmicflows-2 data, in which all galaxies forming a 
group (of two or more) are merged to one data entry. The grouped Cosmicflows-2 
data consists of 4,885 entries. Six methodologies are used for distance estimation: 
Cepheid star pulsations; the luminosity terminus of stars at the tip of the red giant 
branch; surface brightness fluctuations of the ensemble of stars in elliptical galaxies; 
type Ia supernovae; the fundamental plane in luminosity, radius and velocity 
dispersion of elliptical galaxies; and the Tully—Fisher (TF) correlation between the 
luminosities and rotation rates of spiral galaxies. 


Wiener filter and constrained realizations. In the standard model of 

cosmology, the linear velocity field constitutes a Gaussian random vector field”. 
The Cosmicflows-2 dataset, as all other available velocity surveys, provides a 
sparse, incomplete, inhomogeneous and very noisy sampling of the observed flow. 
The Bayesian formalism of the Wiener filter and constrained realizations provides 
the optimal methodology for the reconstruction (estimation) of the underlying 
velocity field and the associated uncertainties in the linear regime'!****. 

The Wiener-filter/constrained-realizations reconstruction is based on an 
assumed prior cosmological model — the ACDM model with cosmological 
parameters inferred from the Wilkinson Microwave Anisotropy Probe (WMAP). 


The current Wiener-filter and constrained-realization fields are the ones reported 
in our bulk velocity article*. The results presented here are insensitive to the exact 
values of the ACDM parameters, in particular to the differences between the 
WMAP and Planck parameters. 


Cosmic V-web. The cosmic web is defined here by the means of the V-web 
model”. The normalized velocity shear tensor at a given grid cell is defined by: 


1 
Zap =— zm oe? + Opa) (1) 


The standard definition of the velocity shear tensor is modified here by the 
Hubble constant (H,) normalization, which makes it dimensionless. The minus 
sign is introduced so that a positive eigenvalue corresponds to a contraction. 
Eigenvalues are ordered by decreasing value, hence ê; points in the direction of 
maximum collapse and ê; points toward maximum expansion. 

The V-web model starts with the continuous velocity field and its associated 
velocity shear tensor (Equation (1)). Consider a given point in space at which the 
shear tensor is evaluated, and thereby its eigenvalues and eigenvectors. The V-web 
is defined by a threshold value (A,,) — a free parameter that defines the web. 

The number of eigenvalues above 4, defines the web classification at that point: 
0, 1, 2 or 3 corresponds to the point being a void, sheet, filament or knot. 

The V-web is defined by the effective resolution of the velocity field and 
by the value of the threshold. Here a Gaussian smoothing of R, = 250 km s™' and 
A» = 0.04 are assumed. 

The sources and sinks of the velocity field, namely the repellers and attractors, 
are closely associated with the voids and knots of the V-web. The voids (knots) are 
regions of diverging (converging) flow, namely regions where the Hessian of the 
velocity potential is negative (positive) definite, yet these regions are in general 
moving with respect to the CMB frame of reference. The repellers and attractors 
are defined here as stationary voids and knots respectively, and hence correspond 
to local extrema of the gravitational potential. Note that the Great Attractor, which 
moves towards the Shapley attractor”, is not an attractor as defined here. 


Multipole expansion of the flow. A first-order expansion of a potential (that is, 
irrotational) velocity field, v(r), around a point labeled by 0 yields: 


Valt) X Voa + (Op Va) TB=Voa— Ho Lap rg (2) 


where v,, and ©, are evaluated at the point 0. This expansion is equivalent to 
a dipole and quadrupole expansion of the (velocity) potential. The flow in a sphere 
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Figure 3 | Aitoff projection in galactic coordinates of the principal 
structures and directions that characterize the flow. These are the dipole 
repeller, the Shapley attractor, the CMB dipole, the bulk velocity and the 
three eigenvectors of the velocity shear tensor. The latter two quantities 
are evaluated across spheres of radii ranging up to 15,000 km s~. The close 
anti-alignment of the bulk velocity with the repeller and the alignment 

of the e, eigenvector with the attractor are robustly manifested. The 
eigenvectors define lines without direction; hence their two endpoints 

are plotted. The anti-apex of the bulk flow and CMB dipole are plotted to 
improve the visualization of their anti-alignment. 


of radius R is modelled here as the sum of a bulk flow, V,,,.(R), and a shear term, 
©,,(R), in the manner of Equation (2). The parameters of the model (that is, 

the bulk velocity vector and the symmetric tensor) are found by minimizing the 
quadratic residual between the model and the actual velocity field, with a 
spherical top-hat window function weighting. 


Streamlines. In the linear regime the flow is irrotational —that is, it is a potential 
flow — and hence the velocity field can be written as a gradient of a scalar 
(velocity) potential, v= V¢,(r). In the linear regime the peculiar velocity (v) and 
ea Q, and 

H are, respectively, the time-dependent cosmological density parameters for 
matter and for dark energy, and Hubble’s constant. The velocity and gravitational 
potential are similarly related. Inspired by the similarity between the gravitational 
potential in linear theory, hence also the velocity potential, to the electrical 
potential in electrostatics, we present the flow field by field lines which we call 
here stream or flow lines. The line equation of a streamline, I(s), where s is the line 
parameter, is calculated by integrating dl(s) = v(I(s))ds. The numerical calculation 
of a streamline involves the determination of the seeds of the streamlines and the 
number of integration steps. For a small number of integration steps and a regular 
grid of seeds, the streamlines resemble velocity arrows. For a large number 

of steps, the flow and anti-flow lines are either trapped by attractors or repellers 
or leave the box. It should be emphasized here that the streamlines are a graphical 
means for the presentation of a vector field and do not represent trajectories 

of objects. Mass elements and galaxies move only a few megaparsecs in a 

Hubble time along streamlines en route from a repeller to an attractor. 


gravitational field (g) are simply related by v= g where Q, 


‘m? 


Quantitative comparison of repeller versus attractor. The contribution 

(in km s~’) of the dipole repeller and the Shapley attractor to the motion of the 
Local Group is compared by means of setting spheres of radius R around both 
objects, and calculating the velocity induced by the density inside these spheres 

at the Local Group. The linear-theory density-velocity relation is used in such a 
calculation. A word of caution is needed before proceeding to the comparison. The 
attractor region is much better sampled by the Cosmicflows-2 data than that of 

the dipole repeller; hence the Wiener-filter suppression of the signal is stronger at 
the dipole repeller, and the constrained variance of the residual around the mean 

is larger there (Supplementary Fig. 1). The figure presents the contribution (in 
kms") to the Local Group velocity projected onto the direction of the CMB dipole. 
Out to Ræ 8,000 km s~’, the contributions of the attractor and repeller are almost 
identical. At 8,000 kms“, the repeller and the attractor contribute 59 + 26 and 

67 +27 km s™', respectively, to the CMB dipole. The slight deficiency of the dipole 
repeller is fully consistent with the enhanced suppression of the Wiener filter. The 
Perseus—Pisces supercluster, at a distance of R > 8,000 km s~! away from the repeller, 
leads to the decrease in the velocity contributed by the dipole-repeller-centric 


spheres. The dip is the reflection of the tug-of-war between the Perseus—Pisces 
supercluster and the Great Attractor’. Very similar results are obtained when 
considering the contribution to the amplitude of the full 3D velocity of the Local 
Group. The conclusion that follows here is that out to R % 8,000 km s~', the Shapley 
attractor’s basin of attraction and the dipole repeller’s basin of repulsion contribute 
equally to the Local Group motion. 


Uncertainties assessment. The probability distribution of the alignment 

of the bulk velocity and the eigenvectors of the shear tensor is sampled by 
means of an ensemble of 20 constrained realizations, constrained by the 
Cosmicflows-2 data and evaluated within the WMAP parameters of the 

ACDM model. The constrained realizations are evaluated on a grid of 

size 256° spanning a box of 256,000 km s™ on its side. The bulk velocity and 

the velocity shear tensor are obtained by a convolution of the velocity field 

with a spherical top-hat window of radius R and are evaluated at the centre 

of the box, in other words the location of the Local Group. Supplementary Fig. 2 
presents the alignment of the bulk velocity and the third eigenvector of the shear 
tensor with the dipole repeller and the Shapley attractor, respectively, over a range 
of radii of R= (20, 30,..., 300) x 100 km s™'. The uncertainty in Hpux(R) is used to 
assess the uncertainty in the position of the dipole repeller. At Rx 6,000 km s7, 
we find /4,,,(R) = —0.96 + 0.04. Assuming the dipole repeller to be responsible 
for the direction of the bulk velocity, the uncertainty in u is translated to a 
projected distance (at a radial distance of 16,000 km s™') of ARpr % 4,500 km s™'. 
The uncertainty in /#,,,, changes significantly over the range R= 11,000 to 
20,000 km s~'. This again translates to an uncertainty in the radial position 

of roughly 4,500 km s™'. 

The robustness of the dipole repeller against subsampling of the data is 
considered. First, the Cosmicflows-2 data is subsampled by type: the ‘singles’ 
subsample consists of 4,264 data entries based on a single galaxy only, and 
the “TF-singles’ subsample is made of 3,943 single data with TF-only distances. 
The ungrouped Cosmicflows-2 data consists of 8,399 galaxies. The relative 
fractional distance error of the TF distances is 0.20, and for all the singles 
is 0.186 + 0.046. The majority of the non-“TF-singles’ are data with supernova 
type-Ia distances for which the relative distance errors are 0.088 + 0.007. 

In assessing the constraining power of the various subsamples, one should consider 
both the number of data points and the magnitude of the errors. For comparison, 
the relative error for the grouped entries (that is, data entries based on more than 
one galaxy) is 0.098 + 0.037. The anti-flow streamlines of the two subsamples 

are shown in Supplementary Fig. 3, and the details of their repellers are given 

in Supplementary Table 1. For both subsamples, a repeller (TF-only) or two 

(all singles) have been found at the back side of the CMB dipole; they both ‘push, 
and the dipole repeller is located within an underdense region. 

Next, a possible ‘edge of the data’ systematic effect is considered. The analysis 
is conducted by trimming the full Cosmicflows-2 dataset in spheres of radii 
6,000, 8,000 and 10,000 km s~}, which contain 49%, 67% and 82% of the full 
data, respectively. The Wiener filter has been applied to these subsets of data, 
and the resulting anti-flow field has been compared with that of the full data 
(Supplementary Fig. 3). The overall structure of the flow fields of the subsamples 
follows that of the full Cosmicflows-2 data. The anti-flow converges into two 
repellers in the 6,000 kms“ case and into single repellers for the 8,000 and 
10,000 km s™ cases. Supplementary Table 1 provides the location and distances 
from the Local Group and the dipole repeller for each case. It also provides the 
fraction of data contained in each subsample. It is remarkable that with less than 
half the data (6,000 kms“ case), the Wiener filter recovers the general anti-flow 
towards the general position of the dipole repeller. Taking the mean position of the 
two repellers, it is found to be a mere 2,500 kms“ from the dipole repeller. Single 
repellers are found for the larger subsamples that consist of 67% and 82% of the 
data at distances of 4,300 and 500 km s™' from the dipole repeller, respectively. 

We conclude that the dipole repeller is not a fictitious structure induced by 
an ‘edge of the data’ effect, and that subsets of the data, chosen either by distance 
or galaxy type, uncover a basin of repulsion that ‘pushes’ the Local Group in the 
direction pointed by the CMB dipole. 


Data availability. The data that support the plots within this paper and other 
findings of this study are available from the corresponding author upon reasonable 
request. They are also available freely at the Extragalactic Distance Database 
(http://edd.ifa.hawaii.edu/) and through the NED interface. Use of the data and 
flow model must cite this article. 
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